Microwave absorption study of pinning regimes in Ba(Fei_xCox)2As2 single crystals 
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Magnetic field dependent modulated microwave absorption (MMWA) measurements have been 
carried out to investigate vortex pinning effects in single crystals of the iron-based high-T^ supercon- 
ductor Ba(Fei-xCox)2As2 with three different cobalt doping levels of a: = 0.07, 0.09, and 0.11. The 
dependence of the MMWA hysteresis loops on temperature, magnetic field, and Co-concentration 
have been measured and analyzed using a theoretical model of microwave absorption in supercon- 
ductors. The analysis reveals that in an underdoped crystal (x = 0.07) the so called 5Tc-pinning due 
to magnetically ordered regions defines the temperature dependence of the critical current density, 
while in the optimally doped {x = 0.09) and overdoped [x = 0.11) samples the pinning is governed 
by structural imperfections due to inhomogeneous distribution of cobalt dopant and has the so called 
51 character. 



I. INTRODUCTION 

The Ba(Fei_xCox)2As2 compound is likely the most 
investigated material among all known superconducting 
pnictides. Such an interest can be explained by both 
unusual physical properties inherent in this compound 
and availability of high quality single crystals. For the 
crystal growth most commonly a self-flux or the Bridg- 
man techniques are usedi. Crystals grown by different 
techniques may vary in properties, in particular, in the 
pinning strength and, as a result, in the superconduct- 
ing critical current density. The vortex pinning effects 
in the Ba(Fei_xCox)2As2 single crystals have been stud- 
ied by transport, magnetic, magneto-optic and relaxation 
measurements (see, for example, Refsi^"— ). The depen- 
dence of the pinning strength on environmental condi- 
tions (temperature T and magnetic field noHa ) and the 
cobalt doping level has been investigated. It has been 
established that the Ba(Fei_xCox)2As2 samples have a 
large critical current density and a steep slope of 
the irreversibility line Hirr{T)^. A symmetric shape of 
the pinning strength function Fp versus field H allowed 
the authors of Ref4 to suggest the presence of a strong 
pinning nanostructure. By simultaneous analysis of the 
magnetization curves M{H) and the domain structure, 
the authors of Ref.^ concluded that the twin boundaries 
play a role of such nanostructure. 

In the present work we have investigated 
Ba(Fei_xCox)2As2 crystals with different Co con- 
centration (from underdoped to overdoped) using the 
magnetic field dependent modulated microwave ab- 
sorption (MMWA) measurements. We have studied 
the evolution of the MMWA signal as a function of 
temperature, applied magnetic field and doping level. 
The relevant superconducting parameters, in particular, 
the critical current density have been estimated through 
numerical simulation of the experimental MMWA loops. 
Furthermore, the MMWA method used in the present 



work has enabled, beyond the scope of the static magne- 
tization measurements, not only to estimate the critical 
current density, but also to obtain additional insights 
into other effects, in particular on the influence of the 
thermal fluctuations and on the pinning center density. 
The theoretical analysis enables a distinction of pinning 
regimes in the studied single crystals with different 
Co-doping, as well as to establish the origin of pinning 
centers. In particular, we have found clear indications 
for pinning of vortices on magnetically ordered nanoscale 
regions in the underdoped sample, whereas at higher Co 
doping the pinning is related to local structural disorder 
and imperfections. 



II. EXPERIMENTAL 

Ba(Fei_xCox)2As2 single crystals with the Co con- 
tent X — 0.07, 0.09 and 0.11 were grown by the verti- 
cal Bridgman technique. The preparation procedure is 
described in detail in Ref.^. The Co concentration was 
determined by elemental analysis using energy dispersive 
X-ray (EDX) mode of a scanning electron microscope 
(see Ref.^ for details). All x values referred to in this 
work are the Co contents as revealed by EDX. Among all 
crystals the one with x = 0.09 has the maximal super- 
conducting critical temperature Tc = 26 K (Fig. 1 and 
Table 1) and thus is considered as an optimally doped 
(OP) sample. Respectively, the crystal with x = 0.07 
is denoted as underdoped (UD), and the crystal with 
X = 0.11 as overdoped (OD). In literature the optimal 
cobalt concentration x for superconductivity in this series 
is reported to lie in the range from 6% to 10% depend- 
ing on the method and regimes of preparation^^"— . Note 
that for Co concentration less than optimal two areas in 
the phase diagram T vs x can be distinguished: the one 
corresponding to the spin density wave (SDW) state (at 
high temperature) and the other corresponding to the su- 
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FIG. 1: Temperature dependence of AC-susceptibility for 
three Ba(Fei_xCox)2 As2 samples with different Co concen- 
tration. 



TABLE I: Superconducting transition parameters for 
Ba(Fei_xCox)2As2 single crystals. 



Sample 


Doping level, x 


r°", K 


rr, K 


AT,, K 


UD 


0.07 


26 


23.8 


5.2 


OP 


0.09 


26.9 


25.9 


1.4 


OD 


0.11 


23.7 


21.5 


3.5 



perconducting state (at low temperatures)^. According 
to a neutron scattering study^, at these concentrations 
at T < Tc the superconducting phase coexists with the 
magnetically ordered phase. 

The parameters of the superconducting transition were 
determined from the temperature dependence of AC sus- 
ceptibility (Fig. 1). The obtained values are summarized 
in Table 1. The onset temperature T°" was determined 
as a temperature at which the AC susceptibility signal 
starts to deviate from the zero line. T™ is the temper- 
ature at the half-height of the transition. The width of 
the transition Tc is determined as the interval between 
the points of 10% and 90% of the total transition height. 
Notably the UD sample has the broadest transition. 

There is no indication of the magnetic phase transition 
in the temperature dependence of the AC susceptibility 
in the studied samples. However for the UD single crys- 
tal it is observed in the resistivity p{T) measurements 
(Fig. 2) which turn out to be a very sensitive probe 
of the magnetic phase transition manifesting in a clear 
minimum of the derivative curve d[p{T)]/dT [for details 
scc^). In fact, following the procedure discussed in Ref,^, 
besides a transition to the magnetically ordered phase of 
SDW type at Tsdw = (40 ± 2) K corresponding to a 
minimum of d[p(T)]/dT we can identify for the UD sam- 
ple a structural transition at Tg — (62 ± 10) K from the 
minimum of the p{T) curve, and finally a sharp super- 
conducting transition onset at Tc = (25 ± 0.2) K (Fig. 
2). 

The MMWA measurements on these samples were car- 
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FIG. 2: Temperature dependence of resistivity (closed circles) 
and its derivative dp/dT (open circles) for the UD sample, re- 
vealing the structural Ts, magnetic Tsdw and superconduct- 
ing Tc transitions. 



ried out in the temperature range of 10 300 K and the 
magnetic field range of 0.75 T. The experiments were 
performed on a Bruker BER-418s ESR-spectrometer at 
a frequency of 9.4 GHz (X-band). The DC field was 
modulated with a frequency of 100 kHz with the ampli- 
tude of 0.1 to 10 Oe. The MMWA signal is detected 
and amplified by the lock-in amplifier at the first har- 
monic of the modulation. The crystal afc-plane was ori- 
ented perpendicular to the applied magnetic field. The 
MMWA signal was registered at fixed temperature upon 
sweeping the magnetic field up and down. We have ob- 
served the field dependent hysteresis of the MMWA sig- 
nal which refiects the magnetic irreversibility due to vor- 
tex pinning (see Fig. 3). The hysteresis amplitude L was 
determined from experimental MMWA loops as the dif- 
ference between magnitudes of the power absorbed upon 
sweeping the magnetic field up {Pup) and down {Pdown), 

T — Pup Pdown • 

The measurement protocol was as follows: A sample 
was cooled from T > Tc down to the measurement tem- 
perature in a fixed magnetic field p^oHres — 5 mT (resid- 
ual field of the spectrometer magnet). Then the applied 
DC field was swept with a rate of 4 mT/s from Hres to 
Umax and back. poHmax was varied from 0.5 to 0.75 T 
depending on a particular measurement. Further de- 
tails on the measurement technique, procedure and the 
physical basis of the MMWA method can be found in 
Refsi^iii. The hysteresis amplitude L appears to be a 
function of cobalt concentration, magnetic field, and tem- 
perature. 



III. RESULTS AND DISCUSSION 

In the temperature range near to the critical temper- 
ature the hysteresis loop closes at a certain field value 
Hirr called the irreversibility field. The temperature de- 
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pendence of the irreversibility field Hirr(T) defines the 
irreversibility line, which separates the areas with zero 
and nonzero values of the superconducting critical cur- 
rent density on the temperature/magnetic field phase di- 
agram. The plot of Hirr versus the reduced tempera- 
ture T/Tc (Fig. 4) allows correct comparison of the irre- 
versibility lines of the samples with different critical tem- 
perature. One can see that the irreversibility line of the 
sample with doping level x = 0.07 (UD) is the steepest. 

The temperature dependence of the amplitude of the 
MMWA hysteresis loop L{T) was obtained at fixed val- 
ues of the applied magnetic field. The L{T) depen- 
dence for three Ba(Fei_xCox)2As2 single crystals is pre- 
sented in Fig. 5 for two selected fields fio^a = O.f and 
0.45 T. These fields are indicated by dotted vertical lines 
in Fig. 3. The field 0.1 T was chosen to be somewhat 
higher than the position of the low-field peak occurring 
at /xpTJa < 0.05 T. This peak arises due to two competing 
effects. The increase of the MMWA signal upon starting 
the sweep of the applied field is due to the rearrangement 
of the vortex distribution. When the field begins to rise, 
the superconducting sample is gradually filled by vor- 
tices, entering the sample from edges towards the center. 
In this regime the growth of the MMWA signal is due 
to the increase of the number of vortices [the term in 
Eq.(l) in fact describes this increase, see below]. As soon 
as the first vortex reaches the center of the sample, the 
volume through which the critical current flows does not 
change anymore. The magnitude of the MMWA signal is 
now determined by the critical current density which de- 
creases with increasing the field strength. Therefore the 
MMWA signal decreases. These two effects are respon- 
sible for a characteristic peak of the MMWA response 
at low fields (Fig. 3). The second (higher) magnetic 
field value of 0.45 T was chosen such, so that to collect 
a sufficient number of experimental data points in the 
irreversible area where L 7^ 0. 

For the UD sample the L first increases with decreas- 
ing temperature, reaches a maximum at T = 21.5 K and 
then decreases. The position of this maximum is indepen- 
dent of the magnetic field value. The L{T) dependence 
for the OD sample is different: Near the superconduct- 
ing transition the amplitude of the loop increases quite 
steeply and tends to level off in the lower temperature 
range. The OP sample reveals a behavior intermediate 
between UD and OD. At /^o-ffa = 0.1 T the L{T) depen- 
dence is similar to that of the UD sample, whereas at a 
higher field /io^a ~ 0.45 T it is similar to that of the OD 
sample. 

The analytical form of L{T) is a complex function of 
the temperature dependent critical current jc , the am- 
plitude u of thermal fluctuations of the flux-line positions 
and the viscosity of the vortex matter 
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FIG. 3: MMWA hysteresis loops for the UD 
Ba(Fei_xCox)2As2 single crystal with x — 0.07 at tem- 
peratures 18 K, 22 K, and 24 K (from top to bottom). The 
loops are shifted on the vertical scale relative to each other 
for convenience. Horizontal arrows show the direction of the 
magnetic field sweep (first up-field followed by down-field). 
The vertical arrow indicates the irreversibility field at which 
the hysteresis loop closes. The vertical dotted lines indicate 
the fields at which the amplitude of the loop L was measured 
as a function of temperature. A small feature at ~ 0.175 T 
corresponds to an EPR signal from iron ions present in the 
glass ampoule in which the sample was placed. 
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FIG. 4: Irreversibility field versus reduced temperature for 
three samples of Ba(Fei_xCox)2As2. 



Here Co is a constant defined by the sample size and by 
the registration conditions of the hysteresis loop, and a 
is the mean distance between pinning centers. Function 
(1) includes integrals over the instantaneous distributions 
of field and current in a sample for sweeping the field 
up {F^{H,j)) and down {F^{H,j)). The distributions 
depend on the shape and sizes of the sample, as well as 
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on the critical current density. A detailed description 
of the i(r, H) function and its components is given in 
Ref. 10, where a theoretical model of the hysteresis of 
the microwave absorption has been developed. 

The 77 (r) dependence can be obtained from the re- 
lation 77 = ^oBc2/Pn : where $0 is the magnetic flux 
quantum, pn is the normal-state resistivity^^. The tem- 
perature dependence of the upper critical field Bc2{T) 
can be derived from the Ginsburg-Landau theory. Then 
one gets 



V = Vo 



1 - (T/T,)^ 
1 + (T/T.y 



(2) 



with 770 = 77(T ~ 0). 

Eq. (1) depends exponentially on the ratio of the mean 
square of thermal fluctuations (u^) to the square of the 
mean distance between pinning centers . According to 
Ref the fluctuation amplitude depends on temperature 
as follows from the equation (m^) = 47rrA^/$Q^^B^/^. 
The general formula of the temperature dependence of 
the field penetration depth has the form A(T) = Ao(l — 
(r/rc)")-i/2. For Ba(Fei_xCox)2As2 the index 771. varies 
from 2 for the UD samples to ^ 2.5 for the OP and OD 
compoundsi^"— . Thus, the temperature dependence of 
the ratio can be written as 



T 



1 - {T/T,Y 



(3) 



Here c' is the parameter that depends on the material 
properties, in particular on the form and on the number 
of the structural defects and impurities, which serve as 
pinning centers. It also depends on the magnetic field, 
specifically at low fields usually c' cx B~^l'^ . 

There is a number of factors that affect the jc(2^) func- 
tion, such as the chemical composition, the material state 
(e.g., poly- or single crystals, or thin films), type and 
concentration of defects, magnetic field strength etcJi. 
It was found experimentally^ that in Ba(Fei_xCox)2As2 
single crystals the critical current density is a power func- 
tion of temperature which can be written as: 



jc = jco(i-(T/r,fr 



(4) 



Here jco is the critical current density at zero temper- 
ature, and "p and n are indices defined by the type of 
pinning centers. This form of the power law extensively 
discussed in RefsJ^"— is characteristic for the most of 
high-Tc superconductors including cuprates. According 
to theoretical studieai^ii^, in low magnetic field these in- 
dexes amount to "p — 2 and n — b/2 for the (5Z-pinning 
inherent in superconductors with weak pinning centers 
of a point type. (The (5^-pinning arises due to spatial 
variations in the charge carrier mean free path I and lo- 
cal lowering of the superconducting order parameter near 
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FIG. 5: Temperature dependences of the MMWA hysteresis 
loop amplitude for the UD (a), OP (b), and OD single crys- 
tals (c) at a magnetic field of 0.1 T (open circles) and 0.45 
T (closed circles), respectively. Solid and dashed curves are 
calculated with Eq. (1). 



lattice defectai^.) The presence of extended defects (in 
particular, twin domain boundaries or impurity phase in- 
clusions) results in the so-called (5Tc-pinning with p = 2 
and n = 7/6. (5Tc-pinning is associated with disorder at 
critical temperature Tc .) 

We simulated the temperature dependence of the hys- 
teresis amplitude using Eq.(l) and taking the tempera- 
ture functions 77(T), {u^{T)) and jc{T) given by Eqs.(2-4) 
into account. The only fitting parameters were jco, c', n 

7 /2 

and Cq = Co/rji^ . The index p was fixed at p = 2 
in accordance with the theory requirementsi^,. We used 
m = 2 for the UD sample and 777 = 2.5 for description 
of the data for OP and OD samples in accordance with 
Refs.— . The results of the simulation are shown in 
Fig. 5 by solid and dashed lines for the applied field of 
0.45 T and 0.1 T, respectively. The obtained fitting pa- 
rameters are presented in Table 2. (The parameter Cq is 
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TABLE 11: Parameters jco, n and c' obtained by fitting Eq.(l) 
to the experimental dependence L{T) at two magnitudes of 
the applied magnetic field noHa ■ 



Sample 


HoHa, T 


jco, A/cm^ 


n 


c', 1/K 


UD 


0.1 


1.7- 10* 


1.3 ±0.2 


(4.5 ±0.5) • 10"* 


UD 


0.45 


1.5 ■ 10* 


1.1 ±0.2 


(5.8 ±0.7) • 10"* 


OP 


0.1 


0.95 • 10* 


1.2 ±0.2 


(6.8 ±0.8) • 10"* 


OP 


0.45 


0.88 • 10* 


2.6 ±0.4 


(4.6 ± 1.3) • 10"* 


OD 


0.1 


1.4 • 10* 


2.9 ±0.5 


(5.5 ±1.2) -10"* 


OD 


0.45 


1.9 • 10* 


1.9 ±0.5 


(9.8 ±2.4) • 10"* 



not given there as it consists of several cofactors which 
cannot be separated.) 

The jco parameter defines the critical current density. 
First of all, it should be noted that its magnitude is in 
a good agreement with the data obtained from trans- 
port and magnetic measurementa^ii. For the UD sample 
it has the temperature dependence described by Eq.(4) 
with the n index close to 1.2 which corresponds to the 
(5Tc-pinning. According to theoretical predictionsi^iiE 
and experimental finding o^^'^" , the (5Tc-pinning occurs in 
a superconductor with large-size defects. In particular, 
such jc{T) dependence is observed in YBCO films con- 
taining large-size (~ 30 nm) inclusions of the impurity 
phase, which provide effective vortex pinning^'^. Thus, 
taking into account that in the UD crystal the value of 
n is close to that of the STc-pmiung regime, a large jco 
value and its weak field dependence can be related to the 
occurrence of large-size pinning centers in this sample. 
Considering the phase diagram of the Ba(Fei_xCox)2As2 
compoundii^ where the coexistence of the superconduct- 
ing and magnetically ordered phases at 2; < 0.8 has been 
observed-^iSi, one can associate such large-size pinning 
centers with regions of magnetically ordered phase which 
should apparently be present in the UD single crystal 
below TsDW = 40 K. In favor of this assumption, the 
study of a similar iron pnictide Bai_xKxFe2As2 by neu- 
tron scattering, muon-spin rotation and magnetic force 
microscopy revealed magnetic (nonsuperconducting) in- 
clusions of the 65 nm size^^. We conjecture that magnet- 
ically ordered domains in our UD crystal have the size of 
the same order of magnitude. 

Another origin of the large-scale defects and their in- 
fluence on the critical current density was discussed in 
Refi^. Twin boundaries were considered as such extended 
defects in Ba(Fei_xCox)2As2 crystals with x < 0.058 
grown by self-flux technique. However, in crystals with 
higher Co concentration no twin boundary structure was 
revealed. For our UD sample with x = 0.07 we therefore 
do not expect twin boundaries to occur. 

The decrease of jco and the change of the n index to- 
wards the value of 2.5 corresponding to the (5Z-pinning 
(at least in the high-field range) in the OP sample (see 
Table 2) is indicative of the pinning due to weak point 
defects^. As suggested in Ref.^^ such pinning centers 



can be formed as a result of inhomogeneous distribution 
of the dopant atoms. It is therefore likely that also in 
our OP and OD samples the Co impurity atoms serve 
as point pinning centers. This assumption is supported 
by the observation of the increase of the critical current 
jco and of the c' parameter, related to the center den- 
sity, upon transition from OP to OD conditions. This 
can be naturally related with the increasing number of 
pinning centers. Some enhancement of the jco value of 
the OD crystal with the applied magnetic field is proba- 
bly connected with the fish-tail effect observed frequently 
in Ba(Fei_xCox)2As2 (see, for example, Refsi^i^). We 
note however that the steepest irreversibility line is ob- 
served for the UD Ba(Fei_xCox)2As2 single crystal (Fig. 
4) suggesting that magnetically ordered phase inclusions 
are most effective pinning centers for superconducting 
vortices. 

Considering the data in Table 2 one should note that 
the parameter c' does not obey the B~^^^ dependence, 
because it is valid in the low field regime only. The c' even 
increases with B for the UD and OD samples. In the lat- 
ter case the increase is most pronounced. For this sample 
also the jc increases which may be related to the fish-tail 
effect, as discussed above. The functional dependence 

< >oc B"^/^ is a consequence of the dependence of 

< > on the shear modulus cge , which in turn is propor- 
tional to ijHiil. This relationship holds well in the low 
field regime. However, with increasing the field two addi- 
tional effects may become relevant. The first effect is the 
softening of the shear modulus which yields a decrease of 
C66 , as discussed in Ref 1^. The second effect pointed out 
in Ref.^^ is the formation of vortex bundles. Both eff ects 
can result in the increase of the vortex shift \/ < > un- 
der the action of thermal fiuctuations, as compared with 
the rigid (elastic) vortex lattice. Plausibly these effects 
may be responsible for an increase of c' and also of jc 
with magnetic field in the case of the OD sample. 

It is possible to estimate t he mean value of the thermal 
fluctuation amplitude u — V < u'^ > and the distance be- 
tween pinning centers a on the basis of the obtained data 
and the known functional dependences. Assuming that 
the above discussed < u^{B) > dependence is valid at 
low fields and taking the value Aq = 2080 we ob- 
tain the following estimates of u and a for the UD, OP 
and OD samples at T = 15 K and H = Q.l T, respec- 
tively: u = 27.7 A, 25 A, 28 A; and a = 263 A, 214 A, 
236 A. One can see that the maximum distance between 
pinning centers, which corresponds to the minimum of 
their density, is revealed for the UD sample. This is in 
agreement with the minimum doping level in this sample 
and supports an assumption on the droplet form of pin- 
ning centers. A somewhat surprising increase of a in the 
OD sample, as compared with the OP sample, suggests a 
possible clustering of dopant atoms. In this scenario the 
concentration of clusters is smaller than that of point de- 
fects. However they are more efficient as pinning centers 
yielding an increase of the critical current density (see 
Table 2). 
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IV. CONCLUSIONS 

A set of the Ba(Fei_xCox)2As2 single crystals with Co 
concentration x = 0.07, 0.09, and 0.11 was studied by the 
MMWA method. The analysis of the temperature and 
field dependences of the MMWA hysteresis amplitude en- 
ables to estimate the critical current density and the crit- 
ical indices of the 3c{T) function, and to determine their 
dependence on temperature and Co doping. On the basis 
of our data we conclude that in the underdoped crystal 
[x = 0.07) the pinning is of a 6Tc type related to the 
presence of regions of a nonsuperconducting (magneti- 
cally ordered) phase which provides the steepest slope of 
the irreversibility line for the underdoped sample. On 
the other hand, in the optimally doped and overdoped 
single crystals the pinning is caused by structural imper- 
fections due to an inhomogeneous distribution of cobalt 



dopant. Because of this, the ^/-pinning is inherent in 
these crystals. 
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